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Abstract 

An observation of autocorrelation of Wilson loops on lattice is presented, espe- 
cially for the small separation in Markov chain. We give a possible explanation 
for such behavior. We also present the dependence of autocorrelation behavior 
on the chosen operators in this paper. 
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1 Introduction 



In order to extract a meaningful physical quantity, e.g. glueball mass, from MC 
simulation of lattice gauge theory, we have to construct the correlation function 
G(t) for chosen operator o(x,t): 

G(f) = £-<o(M)o(0,0) >, 
to extract the glueball mass m when t — > oo, 

G{t) oc e~ mt . 

It is the average over the independent configurations. Meanwhile, we also have to 
make correctly error estimate for the measurement. We should generate a sequence 
of configurations and then choose samples separated by enough sweeps to ensure 
their independence. Therefore it is needed to investigate the correction from auto- 
correlation between configurations in Markov chain. Denoting sequence in Markov 
chain by Tmc, some authorsEJ have studied autocorrelation with different algo- 
rithm and regard normalized autocorrelation asymptotic behavior with separation 
ATmc i n Markov chain as e~ ATMC ^ Texp . In order to improve the fitting, ref. 
redefined r exp via 

AT; 



, . MC 

T exp = iimSUp — : — — r, 

atmc^oo -logp(AT MC ) 

where 

; Am \ < ViVi+ATMc > - < V i >< V i+A T MC > 
Pi^-J-MC) — 

< ViV{ > — < V i >< V i > 

is normalized autocorrelation and v , is the measured value of certain operator in 
the i th configuration. But this formula always works at large ATmc region^. 

Since autocorrelation contribution to error is mainly from small ATmc region, 
we investigates autocorrelation at small ATmc using MC simulation in this paper. 
From the simulation, we find that autocorrelation behavior is like p oc e _Cv/ATMC 
rather than e~ cATMC , especially at small ATmc region. This discrepancy is impor- 
tant since the main autocorrelation contribution to error comes from small ATmc 
region. In section 2, we give a possible explanation for such behavior. 

We also investigate the dependence of autocorrelation with the operators chosen 
for extracting the glueball mass. It is found that autocorrelation depends on the 
chosen operators and we also give a possible explanation of this dependence in 
section 2. We present our simulation results in section 3 and the conclusion in 
section 4. 

We adopt improved action on anisotropic lattice ( one can find the notations 
and details in ref. B): 



H ~ Pi 3^? s l^T 12£«2 12^p' 
To simplify, we restrict ourselves to make the simulation in SU(2) pure gauge fields. 
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2 An Explanation for the Autocorrelation Be- 
havior at Small ATmc Region 

A configuration can be described by AN 4 SU(2) matrices, where A 4 is lattice site 

number. This means we can depict configuration by 12iV 4 = M variables, a factor 
of 3 multiplied here is due to 3-dimension Lie algebra for each link(We regard it 
as E 3 space for simplification). Therefore, configuration can be regarded as one 
point P in M- dimensional Euclidean space S. To simplify, we assume the spacing 
of time-like and space-like links are of the same. Then all the M parameters, or 
coordinates, are of equivalence due to the periodic boundary condition. 

When updating configurations, point P forms a trajectory in S space. Define 
distance between two point x — (xi, x 2 , ■ • ■ , xm),U = (yi, 2/2, • • • , Vm) £ S by 



In one sample, we regard the same-time-slice sub-configuration as a point Pi in 
9A 3 (= Mi ) -dimensional space S\( There are N such sub-configurations in one 
sample. ). As expected, the series of A points in one sample can be regarded 
as the evolution of the sub-configuration with Euclidean time t. All the links in 
these sub-configurations are space-like, and we can define distance r' of two points 
in Si in the similar way. After doing these, we find that point Pi in Si has two 
independent evolution behaviors: the first one is the evolution with T MC at fixed t 
and the second one is the evolution with t at fixed T MC ;(in the same sample). 

Suppose configurations have approached to equilibrium after enough preup- 
datings. Updating configuration n times, we get representation points {x^}(i = 
1,2, •••,n) in S space. We denote the displacement of two adjoining points by 

— s.(i) 

r = x^ +1 ) — x^\ (xW is also a vector.) When n — > oo, we get a sequence of 
r with certain distribution. Due to ergodicity and periodic condition, it can be 
expected that the distribution of r only depends on its modular value r. Since the 
procedure of updating configuration, or Markov chain, forms a trajectory in space 
S, and the procedure is of random, we may regard the procedure as random walk- 
ing in space S. Then mean distance of two points separated by AT MC in Markov 
chain is given by < r >= y / 2D AT MC , where D is diffusion coefficient which only 
depends on the choice of algorithm. 

Now let's consider sub-configuration evolves with Tmc at fixed t\ ( the sub- 
configuration consists of space-like links in time-ti-slice). We denote corresponding 

d) ^ W 

points and displacements in Si with {x } and r' . When n is large enough, we 
expect that mean distance of two sub-configurations adjoining in Markov chain is 



yj Mi/ 'Mr where r is mean distance of two representation points in S space ad- 
joining in Markov chain. So, at certain ti, mean distance of two sub-configurations 
separated by AT MC in Markov chain is 



r(x,y) = \j (x~i - y x ) 2 + (x 2 — y?) H h (x M ~ Vm) 2 - 





(1) 
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The distribution of r' has no dominant direction. 

Then we consider sub-configuration evolves with time t in one configuration. 
When configuration approaches to the equilibrium, the most possible contribution 
to configuration is from the vicinity near the solution of classical equation of motion, 
and we expect the distance of two sub-configurations separated by time At is a well- 
defined function. Up to the first order, mean distance of two sub-configurations 
separated by At in the same sample in Si space can be expanded as 

< r' >~ c\ + c' 2 At, (2) 

where c' 2 > 0. The distribution of r' has no dominant direction, too. Eq.(2) can 
also be written as 

At = ci + c 2 < r' > . (3) 

Since we are interested in operators for extracting glueball mass, we will only 
consider the operator with certain quantum number J PC . Suppose o is such kind 
of the operator and its measurement o(t) in sub-configuration is a function of point 
Pi in Si. Autocorrelation of o(t) is actually the correlation function G(t). One can 
extract the glueball mass by its asympototic behavious 

p{At) oc e- mAt . (4) 

On the other hand, due to Eq.(3), it is 

p(At) oc e - C2m<r,> . (5) 

When we calculate correlation function with certain AT M c or At, we should average 
it in many samples, or average it in subspace Si due to ergodicity. Therefore, p is 

function of r'. Meanwhile, the distribution of r' has no dominant direction with 
certain ATmc and At, so p is only the function of r'. 

Let us, now, consider the evolution of sub-configuration with ATmc again. Due 
to Eqs. (1), (5), we get 

p{AT MC ) oc e -™*V2S!ZEF ( 6 ) 
In our simulation, we find it is better to express the autocorrelation p by 

P(ATmc) OC e -(^+c 2 mWAT MC _ (7) 

Presumably, it is due to finite size effect in lattice simulation. One can find that 
autocorrelation gets weaker with the increasing of the mass extracted from the 
correlation function of the operator G(t)(c2 > 0) . 
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3 Simulation Results 



As we know, one can construct certain combination of Wilson loops to extract the 
mass of J PC glueball states. These special combination of Wilson loops transforms 
according to certain representation of cubic group(Ai, A 2 , E, T 1; T 2 ) §0. We 
choose some of these operators to observe the autocorrelations. These operatore 
are the combination of the certain prototype Wilson loop which are shown in Fig. 

1. Operator a) belongs to representation Af(J p = + ); b) belongs to representa- 
tion E + (J P = 2 + ); c) belongs to representation T 2 + (J P = 2 + ) and d) belongs to 
representation Ai(J p = 0~). The prototype of operators a) and b) are plequettes. 
We choose different combination coefficent and different orientation to make op- 
erator a) belongs to representation Af and operator b) belongs to representation 
E + (J P = 2 + ). Therefore, one can extract different mass states from them. Opera- 
tor b) and c) belong to different representations of cubic group but correspond to 
the same continuum J p states ( Masses extracted from operator b) and c) approach 
to the same value in the continuum limit). 

Our simulation is performed on a 8 3 x 24 anisotropic lattice with £ = 2a = 
3.0 and (3 = 1.0. The mass ma t extract from Af , E + , T 2 + , A~{ operators are 
0.74(2), 1.27(10), 1.21(7), and 1.59(7) respectively. 

We present the operators autocorrelation behavior and their fitting curve in Fig. 

2. Except at the large ATmc region where the noise will overwhelm the signal, the 
formula (7) works very well. 

When analyzing autocorrelation data for operator Pss, spatial plaquette and 
for operator Rss, spatial 2x1 rectangle from ref. 0, one also find our formula (7) 
does work well. Our argument is also supported by the data for different operators 
of ref. ( their Fig. 2), although their data is to describe fermions. 

When ATmc is not small, the space S can not be regarded as Euclidean one, so 
Eq.(l) is not valid. At the same time, the expansion in Eq.(2) is doubtful at large 
At. Therefore, the explanation in section 2 is not advisable at large ATmc region 
or weak autocorrelation region. One may expect that when ATmc is very large, the 
autocorrelation behavior will transform into the behavior e~ ATMC ^ Texp . This can be 
seen obviously in Fig. 2 a). (To guide our eyes, we extract it in Fig. 3 with — Inp 
verse AT M c )• We find that the fitting curve is —Inp = —0.0786 + 0.675V AT M c ( 
for small AT M c) an d —Inp = —2.01 + 0.261ATmc (for large AT M c) respectively. 
But this transition does not happen suddenly. In fact, at intermediate ATmc 
region, the formula — Inp = —1.0 + 0.104ATjy/c works as well as the first formula. 

The dependence of autocorrelation on the mass mat extracted from the corre- 
lation function G(t) is very intricacy. Fitting autocorrelation by p oc e _Cv ' ATjv - fC ', 
Fig. 4 presents c verse ma t by c = —0.65 + 1.94ma t . Our data implies that with the 
increasing of the mass, the autocorrelation gets more weak, as expected in eq.(7). 
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4 Conclusion 



From our simulation, one can find that autocorrelation behavior is like e~ Cv/ATMC 
at small AT MC region. As AT MC increases, autocorrelation approaches to e~ cATMC 
gradually, as one expected. Therefore, we can get fair error estimation from the 
formula e~ Cy/ATMC . An explanation is also given in this paper for such behavior. 

Our simulation also show, the autocorrelation gets more weak with the increas- 
ing of the extracted mass mat- It is consistent with the fact, as one usually expect, 
the autocorrelation is dependent on the operator chosen. Since ++ glueball is the 
lightest state of glueballs, therefore, the correction from autocorrelation is larger 
in ++ glueball case then the others. 
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Figure Captions 



• Figure 1 Prototypes of the corresponding operators a) - d) in the paper. 
Each operator is the different combination of the prototype 00 with direc- 
tions in 3- dimensional space, a) and b) are both plaquette, but we choose 
different combination coefficients to make operator a) belong to representa- 
tion Af and operator b) belong to E + in the cubic group. We also choose 
proper coefficients to make operator c) belong to T 2 + and operator d) belong 
to Ai in the cubic group. 

• Figure 2 autocorrelation behaviors p of operator a) - d) in Fig. f and their 
fitting curves. Fitting curves (dashing curve): a) is p = e - 0786 - - 675 v / AT M c • ^ 

is p = e°-92-2.0 2v /2TW. c ) is p = e 0.7-1.837STW^. and d ) is p = ^.64-2.27^27^ 

The solid curves in a) are p = e 2 01 -°- 261AT McQ ar g e ATmc region) and p = 
e -i.o-o.io4AT MC ^ n ^ ermec ^ a |- e ATmc region) respectively. The vertical coordi- 
nate is p and the horizontal one is ATmc- 

• Figure 3 —Inp verse ma t . The dashing curve is —Inp = —0.0786 + 
0.675 VATa/c). The solid line are -Inp = -2.01 + 0.261AT MC ( large AT MC 
region) and — Inp = 1.0 + 0. 104 ATmc ( intermediate ATmc region) respec- 
tively. The vertical coordinate is —Inp and the horizontal one is the mass 
extracted from the operator. 

• Figure 4 c verse ma t . The fitting line is c = — 0.65 + 1. 9Amat( The dashing 
line is c = 1.487maj). The vertical coordinate is c and the horizontal one is 
the mass extracted from the operator. 
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